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ABSTRACT 

We present XMM-Newton X-ray observations of the interacting galaxy pairs 
NGC 7771/7770 & NGC 2342/2341. In NGC 7771, for the first time we are able to re- 
solve the X-ray emission into a bright central source (Lx ~ 10 41 erg s" 1 ), plus two bright 
(Lx > 10 40 erg s _1 ) ultraluminous X-ray sources (ULXs) located either end of the bar. The 
soft emission of the central source is well-modelled by a two-temperature thermal plasma with 
fcT=0.4/0.7 keV. The hard emission is modelled with a flat absorbed power-law (V ~ 1.7, 
Nh ~ 10 22 cm -2 ), and this together with a low-significance (1.7er) ~ 300eV equivalent 
width emission line at ~6keV are the first indications that NGC 7771 may host a low- 
luminosity AGN. For the bar ULXs, a power-law fit to X-l is improved at the 2.5cr level 
with the addition of a thermal plasma component (kT ~ 0.3 keV), while X-2 is improved 
only at the 1.3cr level with the addition of a disc blackbody component with Tj„ ~ 0.2 keV. 
Both sources are variable on short time-scales implying that their emission is dominated by 
single accreting X-ray binaries (XRBs). The three remaining galaxies, NGC 7770, NGC 2342 
and NGC 2341, have observed X-ray luminosities of 0.2, 1.8 & 0.9 x 10 41 erg s -1 respectively 
(0.3-10keV). Their integrated spectra are also well-modelled by multi-temperature thermal 
plasma components with fcT=0.2-0.7keV, plus power-law continua with slopes of T=1.8- 
2.3 that are likely to represent the integrated emission of populations of XRBs as observed 
in other nearby merger systems. A comparison with other isolated, interacting and merging 
systems shows that all four galaxies follow the established correlations for starburst galaxies 
between X-ray, far-infrared and radio luminosities, demonstrating that their X-ray outputs are 
dominated by their starburst components. 
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1 INTRODUCTION 

It is well known that starburst activity can be triggered, and 
possibly enhanced, by galactic encounters/collisions that may re- 
su lt finally in a mer g ing sy stem. The original N-body studies 
of iToo mre & Toomre 1 19721) and many subsequent studies (e.g. 
Barnes & Hernauist 1996) demonstrated how interactions can re- 
distribute substantial quantities of material towards the central re- 
gions of the galaxies, triggering and fuelling violent bursts of star 
formation. Such encounters also create complex structures such as 
tidal tails, gas bridges and other disrupted morphologies, which 
are easily o bservable in the optical and radio wavelengths (e.g. 
Iffibbard & van Gorkoml~9 96). and can influence the properties of 
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all gas phases from cold neutral gas, to hot ionised gas. Another 
possible consequence of tidal interactions is the formation of a 
galactic bar re sulting from instabilities created in the disc (e.g. 
Noauchi 1988; Barnes &Hernauisll99ll) . Bars are believed to re- 
sult in an inflow of gas to the central regions of the galaxy by the 
extraction of angular momentum from the gas in the disc through 
gravitational torques, thereby fuelling a starburst and/or an active 
galactic nucleus (AGN) iBarnes & Hernquisll99ll) . 



Substantial observational evidence for enhanced star for- 
mation in interacting syst ems has been found a t optical and 
infrare d (IR) wavelengths jLarsorL&_T jnslej Il978t I Jose ph et all 



11984 Lonsdale. Persson. & Matthews! 


1984 iBushousd 
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Kennicutt et al. 


1987; 



ISanders & Mirabel Il996l) . Although enhanced star formation 
activity is found in the discs of interacting galaxies, the ef- 
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feet is much stronger in the nuclear regions. Galaxies that are 
in close pairs or possess a bar show enhanced central radio 
emission, stronger than those of isolated/unbarred systems (e.g. 
Hum mel etaiTl 990). plus high nuclear emission-line luminosities 
and ionization levels I Kennicutt & Kee] 1 1984 iKeel et alj 11985b 
Beravall. Laurikainen. & Aaltol 120031) . It is not currently clear 
whether galaxy-galaxy interactions are sufficient for triggering 
AGN activity, or wheth er the inc reased activity is entirely starburst 
in nature. For example. ISchmittl 12001 ) found that Seyfert galaxies 
have the same percentage of companion galaxies as other activity 
types (e.g. LINERs, HII & transition galaxies), though this could 
simply mean that we observe galaxies going through different 
stages of activity following the movement of gas to the central 
regions, as a delay is expected between the onset of an interaction 
(and an initial starburst phase), and the accretion of the material 
onto a nuclear black hole. Indeed, significant correlations between 
AGN activity and interactions are only found in high mass accre- 
tion rate sy stems such as very luminous/radio-loud quasars (see 
|jogeel2004l for a recent review). 

X-ray studies offer a powerful means to understand the na- 
ture of starburst galaxies (SBGs), because enhanced X-ray emis- 
sion is a manifestation of both star formation and AGN activ- 
ity. With the launch of ROSAT and ASCA, it became possible 
to obser ve the X-ray emission from int eracting/merging galax- 
ies (e .g. |Kbbmno J _Schw^izer J &Mackie] [l997|; |Read_&P^nmad 
ll998tlAlonso-Herrero et alJll999l : lHenriksen & Cousineaull 199*911 
However, the poor spatial resolution and relatively low effec- 
tive area of these missions has hampered attempts to observe 
the X-ray morphology and detailed spectral properties of these 
systems. Now, with the combined spectral and imaging capa- 
bilities of XMM-Newton and Chandra, we have the opportu- 
nity to extend this work to a larger sample of systems at 
different stages of interaction. Many XMM-Newton and Chan- 
dra studies have been made in recent years of merging/closely 
intera cting systems e.g. N GC 4438/39 (The Antenn ae, e.g. 
Fabbiano. Zezas. & Murravl200lh. NGC 4485/4490 jRoberts et all 
20021) , NGC 46 76 (The Mice.lReadl2003T) and post-me rger galaxies 
e.g. NGC 3256 fcira et all2002Ujenkins et all2004bl) and Arp299 
IZezas. Ward. & Murravl2003h . Such studies have revealed that the 
hard (2-10 keV) X-ray emission of SBGs are domina ted by a hand- 
ful of bright extra-nucle ar point sources (see also iKilgard et alJ 
l2002HColbert et aljE o04). the positions of which generally corre- 
late with regions of massive star formation. The brightest of these 
sources are the ultraluminous X-ray sources (ULXs), with X-ray 
luminosities > 10 39 ergs -1 , exceeding the Eddington limit for 
accretion onto a 1.4A/q neutron star by a factor of ~ 10 (see sec- 
tion l5.1.1> . 

Although many merging/post merger galaxies have been stud- 
ied in detail, galaxies in earlier stages of interaction have recently 
been largely neglected. In this paper, we present XMM-Newton X- 
ray observations of NGC 7771/7770 and NGC 2342/2341, two 
pairs of SBGs that appear to be in an early stage of interaction. 
This paper is set out as follows. In section |2| we summarise the 
multiwavelength properties of the two systems. In section [3] we 
outline our X-ray data reduction and analysis methods, and in sec- 
tion|4|we present detailed results of the spectral fitting and timing 
analyses, comparing our results with previous X-ray observations 
of these systems. Finally, we discuss our results in section|5|in the 
context of the properties of other interacting galaxies. Throughout 
this paper, we adopt distances to the galaxies assuming a value of 
J ff =75kms- 1 Mpc -1 . 



2 THE GALAXIES 
2.1 NGC 7771/0 

NGC 7771 is a highly inclined (/~70°) early-type barred spiral 
galaxy (SBa) at a distance of 56.7 Mpc. It is part of the Lyon 
Group of Galaxies No. 483 consisting of NGC 7769, N GC 7771, 
NGC 7786, NGC 7770 and UGC 12828 <Garcidll993l) . and has 
also been shown t o be interacting with a small edge-on disc 
galaxy NGC 7771 A iNordgren et al jl997l) . NGC 7771, NGC 7770 
and NGC 7769 are also catalogued as an interacting triplet 
I Karachent sev et alJll9 88). although in this study we concentrate 
on the closely interacting pair NGC 7771/0, which have an angu- 
lar separation of only ~1.1 arcminutes on the sky corresponding 
to a projected physical separation of 18 kpc at 56.7 Mpc. Note that 
NGC 7769, located ~ 5.4 arcminutes north-west of NGC7771, is 
detected as an extended X-ray source in this XMM observation and 
will be included in an upcoming study. Of the other members of 
the group, NGC 7771 A is not detected, and both NGC 7786 and 
UGC 12828 are outside the field-of-view of the observation. 

NGC 7771 has a high far-IR (FIR) luminosity 
(logL FIR /L Q = 11.1. ISandersetalj Eo03). classifying it as 
a luminous IR bright galaxy (LIRG). Previous studies in the 
optical, radio and IR show that NG C 7771 is going through a 
period of intense star formation ( Smith et alJ I 19991 and references 
therein). The nucleus has be en classified a s an H H-type from 
optical emission line ratios IVeilleux et alJ Il995t) . and shows 
no evidence of broa d Ha emission expected from a n AGN 
IVeilleux et all 1 19951 iDavies. Alonso-Herrero. & Ward Il997h . 
NGC 7771 has a d isturbed morphology and a skewed Ho rotation 
curve <Smithetaljll999h . both of which may be attributable to 
its interaction with N GC 7770. The central s tarburst region has 
been resolved at radio <Neff& Hutching J 1992h and near-IR (NIR) 
ISmith et alJll999l) wavelengths into a nuclear source surrounded 
by a circumnuclear starburst ring ~ 6 arcseconds (1.6 kpc) in 
diameter. The ring is composed of radio-bright and NIR-bright 
clumps, with spatial displacements between the peaks of radio and 
NIR emission indicative of mul tiple gene rations of star formation 
spanning ~4-10Myr ISmith et alJ Il999h . NGC 7771 has been 
previously observed in X-rays using both the Einstein Observatory 
I Fabbi ano. Kim. & Trinchieri| | l992t iRephaeli. Gruber. & Persiol 
1 19951) and ROSAT iDavies et al J 19971 see section l4~T31 . 



2.2 NGC 2342/1 

NGC 2342/23 4 1 is a more distant system at 71 Mpc 
iGoldadereta?] Il997ii) . and is included in the Catalogue 
of Isolated Pairs o f Galaxies in the Northern Hemisphere 
iKarachentsevI fl972l) . The galaxies have an angular separation 
of ~ 2.5 arcminutes, which corresponds to 52 kpc at 71 Mpc. 
This is not a particularly well studied pair, but both galaxi es have 
high IR luminosities of logLFia/ Lq = 10.8 ISanders et alJ 
2003). The larg est galaxy of the pair is NGC 2342, classed as 
S(pec) Ide Vaucouleursl Il99ll) . It shows no evidence of AGN 
activity at any wavelength, and lHo. Filippenko. & Sargent! IT997b) 
classify the nucleus as HII rather than an AGN based on optical 
emission line ratios. This is also suppo rted by the la ck of broad 
Ha emission expected from an AGN lHoetalJll997l). and the 
lack of a strong NIR continuum <GoldaderetalJll997bl) . This 
system was previously observed in the X-ray regime with the 
ROSAT PSP C, and is included in the X-ray survey of galaxy pairs 
by Henriksen & Cousineau 1999 (see section l4.2.4t . 
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Table 1. XMM-Newton observation details. 



Target 


Observation ID 


Date 


Filter 


Useful exposure (s) 


Source count rate (ct 


s -1 , 0.3- 


-lOkeV) 






(yyyy-mm-dd) 




PN MOS 


PN 


MOS 




NGC 7771 


0093190301 


2002-06-21 


Medium 


27148 31343 


7.7xl0~ 2 


2.3x10" 


-2 


X-l 










9.0xl0"" 3 


2.4x10" 


-3 


X-2 










8.1xl0~ 3 


1.8x10" 


-3 


NGC 7770 










2.7xl0~ 2 


7.2x10" 


-3 


NGC 2342 


0093190501 


2002-03-25 


Medium 


24103 25096 


8.7xl0- 2 


2.7x10" 


-2 


NGC 2341 










4.7xl0~ 2 


1.4x10" 


-2 



3 OBSERVATIONS, DATA REDUCTION & ANALYSIS 
TECHNIQUES 

3.1 XMM-Newton Observations & Reduction 

The details of the XMM-Newton observations are shown in TableQ 
NGC 7771/0 was observed with XMM-Newton for ~31 ks on the 
21st June, 2002; NGC 2342/1 for ~29ks on the 25th March, 2002. 
Both observations were performed with the EPIC MOS-1, MOS- 
2 & PN cameras in "Prime Full Window" mode with medium fil- 
ters. The event lists were pipeline-processed using the SAS (Science 
Analysis Software) v5.3.1/5.3.2 for NGC 7771/0 and NGC 2342/1 
respectively, and all data products (spectra, images & lightcurves) 
were created using SAS v5.4.1 using EPIC calibration as of Febru- 
ary 2003 (NGC 7771) and September 2003 (NGC 2342). The data 
were filtered using standard patterns corresponding to single and 
double pixel events for the PN (0 & 1-4), and patterns 0-12 (sin- 
gle to quadruple events) for the MOS, together with f lag=0 to 
remove hot pixels and out-of-field events. 

Both observations were checked for periods of high back- 
ground soft proton flares by accumulating full-field light curves. 
There was no flaring present in the NGC 7771/0 observation, hence 
all —31 ks of MOS-1 & MOS-2 and ~27 ks of PN data were used 
in the subsequent analysis. A small amount of flaring was present 
in the NGC 2342/1 observation, and the times of three peaks (each 
with count rates > 50 counts per second in the 0.3-10 keV PN 
data) were screened from subsequent data analysis, leaving net 
good times of~25 ks (MOS-1 & MOS-2) and ~24ks (PN). 

3.2 XMM-Newton Data Analysis 

We used the SAS task ESPECGET to simultaneously extract source 
and background spectra as well as create response matrices (RMFs) 
and ancillary response files (ARFs) for each source. The spectra for 
NGC 7771 and NGC 7770 were extracted using circular regions 
of 25 arcseconds radius, while the two discrete sources X-l & X- 
2 were extracted using circular regions of 12.5 arcseconds radius 
each. These regions were chosen to include the maximum emis- 
sion from each source while avoiding contamination from nearby 
sources, and all sources were treated as point sources so that point- 
spread function (PSF) encircled energy corrections were applied. In 
each case, the background regions were taken close to each source, 
using as large an area as possible with approximately the same 
DET-Y distance as the source region (to ensure similar low-energy 
noise subtraction). For NGC 2342/1, the spectra were extracted us- 
ing circular regions of 35 and 30 arcseconds radius respectively, 
with background regions selected in the manner just described. The 
resulting spectra were binned to a minimum of 20 counts per bin in 
order to optimise the data for \ 2 statistics. 

The spectral analysis of the XMM data in this study has been 



performed using XSPEC vl 1.1/1 1.3. All errors are given at the 90 
per cent confidence level unless stated otherwise. The MOS and 
PN spectra have been fitted simultaneously in the 0.3-10 keV band 
for each source, and we have included a free normalization constant 
to account for differences in the flux calibration of the three EPIC 
cameras, which differ by < 15 percent. 

The spectral models we use to fit the data of the central 
source in NGC 7771 and the integrated emission from NGC 7770, 
NGC 2342 and NGC 2341 are a power-law continuum (PL) 
to represent the combined emission from compact sources, and 
an optically-thin thermal plasma (MEKAL, also denoted M) to 
model the thermal components expected in SBGs such as galac- 
tic winds and supernova remnants (SNRs). These types of mod- 
els have been used successfully to describe the broad-band X- 
ray emiss ion detected in other SBGs with ASCA (NGC 253 
& M82; IPtaketalJ Il997l). an d XMM-Newton/Chandr a (e.g. 
NGC 253, IPietsch et alj boOll IStrickland et alj l2002t M82, 
I Stevens. Read. & Bravo-GuerrerofcOOl NGC 3256 & NGC 3310. 
iLira et all2002Lljenkins et a l. 2004b). For this study we have cho- 
sen to fix the abundance parameters in the MEKAL models to 
solar values, as previous results have shown that fitting compli- 
cated multi-temperature thermal plasmas with this type of sim- 
plistic spectral model tends t o result i n unrealistically low metal 
abundances of ~0.05-0.3Z© iStrickland & StevensfeOOd and ref- 
erences therein). For the discrete sources in NGC 7771, we have 
also applied models typically used to model emission from X-ray 
binaries (XRBs) such as the multicolour disc blackbody (DISKBB) 
model, which describes thermal emission from an an accretion disc 
jMitsuda et aflll984l : iMakishima et alJll986t) and a conventional 
blackbody (BBODY) model. 

4 RESULTS 
4.1 NGC 7771/0 

4.1.1 X-ray morphology 

FigureQshows soft (0.3-2keV, top) and hard (2-10keV, middle) 
band XMM-Newton EPIC images of the the NGC 7771/0 system. 
The images are the sum of all 3 detectors, and are convolved with 
a lcr (4 arcsecond) 2-D Gaussian mask and displayed with a log- 
arithmic scale. The intensity minima/maxima correspond to sur- 
face brightnesses of 1.2 x 10~ 5 /4.4 x 10~ 4 ct s _1 (soft) and 
8.8 x 10~ 6 /3.5 x 10~ 4 ct s" 1 (hard). The contours correspond 
to 0.4,1.1 & 2.8 x 10~ 4 ct s _1 (soft) and 0.6 & 2.8 x 10" 4 ct s _1 
(hard). Figure[T](bottom) shows a DSS blue image of NGC 777 1/0, 
overlaid with broad-band (0.3-10 keV) X-ray contours correspond- 
ing to surface brightnesses of 0.6,1.4,2.2 & 9.4 x 10~~ 4 ct s _1 . 

The soft X-ray emission is prevalent throughout the central 
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Figure 1. XMM-Newton EPIC soft (0.3-2 keV, top) & hard (2-10 keV, mid- 
dle) X-ray images of NGC 7771/0 with surface brightness contours over- 
laid. Bottom: DSS optical blue image of NGC 7771/0 with broad-band 
(0.3-10 keV) X-ray contours overlaid. The cross marks the optical centre 
of the galaxy. See text for further details. 



and bar regions of NGC 7771. For the first time, we are able to 
resolve the hard X-ray emission into a bright central source plus 
two sources located west and east of the central region. The DSS 
image (FigureQ bottom) shows that the X-ray emission of the cen- 
tral source correlates well with the optical centre of NGC 7771. 
The best available (J2000) position for the centre of NGC 7771 is 
a=23h 51m 24.8s, <5=+20 d 06m 42.6s ±1.2 5 arcseconds from the 
2MASS All-Sky Catalog JCutri et alj|2003l) . and is marked with 
a cross on the OM image. This agrees to the centroid of the X- 
ray emission (a=23h 51m 24.7s, <5=+20d 06m 42.2s) to within the 
XMM positional error of ±1.5 arcseconds (which includes both sta- 
tistical and systematic pointing errors). 

The two off-nuclear sources in NGC 7771 are offi- 
cially designated XMMU J235 122.3+200637 (west) and XMMU 
J235 127.2+200652 (east), although in this paper they are referred 
to as X-l and X-2 respectively. They are located at either end of the 
galactic bar at distances of ~34 arcseconds (X-l) and 36 arcsec- 
onds (X-2) from the nucleus, which places them almost exactly at 
the corotation radius of the bar (35 arcseconds. ISmith et all 1999). 
We have checked whether the central and bar sources are extended 
by fitting their radial profiles in the stacked (PN+MOS) images us- 
ing the IRAF task imexam, and comparing their full-widths at half- 
maximum (FWHM) to that of the XMM EPIC PSF (~ 6 arcsec- 
onds). Both the nuclear source and the bar source X-l show evi- 
dence of extended emission, especially in the soft 0.3-2 keV band 
where the FWHM for each measures ~ 10— 1 1 arcseconds. X-2 may 
also be slightly extended with a FWHM of ~ 8 arcseconds. 

The soft (0.3-2 keV) X-ray emission from the dwarf compan- 
ion NGC 7770 correlates well with the optical image, although it 
does not appear to be a strong source of hard (2-10 keV) emission 
(this agrees with its steep spectral shape, see section l4. 1.3> . There 
is evidence that the soft emission is extended, with a FWHM of up 
to 12 arcseconds in the soft (0.3-2 keV) EPIC images. In addition, 
there appears to be low-level soft diffuse X-ray component extend- 
ing from the south of NGC 7771 towards NGC 7770 that could be a 
consequence of tidal interactions between the galaxies. The lowest 
surface brightness contour that traces this feature in FigureQ(top 
left) represents flux level ~ 5 times that of the average background 
flux level in that area. However, the distance between the peak of 
the emission in both galaxies is ~ 1 arcminute, and we cannot rule 
out that this is an artifact of the large PSF of the EPIC instruments 
(~ 45 arcsecond 90 percent encircled energy radius at 1.5 keV). A 
deep, high-resolution Chandra observation is required to ascertain 
whether this feature is real. 



4.1.2 NGC 7771 Spectral Analysis 

We began the X-ray analysis of the central source in NGC 7771 
by fitting the spectra with single- and two-component spectral 
models with combinations of power-law and MEKAL components, 
each with separate absorbing columns. All single-component mod- 
els were rejected due to unacceptably high values of chi-squared 
(xi > 2), but good fits were obtained with two-component mod- 
els. The soft component can be represented by a warm thermal 
plasma (kT ~ 0.5 keV), and the hard component is well mod- 
elled by either a hard power-law (T ~ 1.5) or a hot thermal plasma 
(kT ~ 19 keV). Although these fits are statistically indistinguish- 
able with xi ~ 0.98 in all cases, the most realistic interpreta- 
tion of the data is the M+PL model, as we know that the hard 2- 
lOkeV emission in nearby SBGs is dominated by a small number 
of bright discrete sources which are likely to be XRBs. However, 
even though hotter gas (typically a few keV) from SNRs is also de- 
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Table 2. Spectral fitting results for NGC 7771/7770. 





PL 


r 


Mi M 2 
N H a kT(keV) N H a kT{keV) 


X 2 /dof 


Fx" 


Lx c 
Obs Unabs 


NGC 7771 Central Source 



Model 1: M+PL (WABS*(PO+MEKAL)) 
, 07 +0.20 1 45+0 06 . 05 , 

lu/ -0.21 lqo -0.08 ' u;l1 

Model 2: M+PL (WABS*PO+WABS*MEKAL) 



-0.06 
-0.05 



1.3912;™ 149 



+0.12 
0.13 



<1.01 



0.53 



+0.08 



Model 3: M+M+PL (WABS*PO+WABS*MEKAL+WABS*MEKAL) 



10.17 4 



1 «i£+0.13 
L66 -0.28 



<0.43 



0.35^ 



6.57 



+2.50 
2.24 



"•"'-0.09 



152.3/156 2.81 



+0.18 

-0.25 



151.5/155 2.80 



+0.11 

0.34 



141.1/152 2.78 4 



1.08 



+0.07 
0.10 



+0.04 



1.08 



'•"'-0.33 



1.28 1 



1.25 



2.24 4 



+0.02 



NGC 7770 



Model 1: M+PL (WABS*(PO+MEKAL)) 

24 ,+0.42 
z "^'-0.32 



Q72 +0.43 



t 



0.50 



+0.07 
0.14 



Model 2: M+PL (WABS*PO+WABS*MEKAL) 
s+0.23 o.n+2.10 



< 1.18 



2.16" 



0.39 



-,+0.13 



Model 3: M+M+PL (WABS*(PO+MEKAL+MEKAL)) 
U6lJ;g 2.26l°;« f 0.15(<0.23) 



52.3/55 



49.7/54 



0.59 



+0.06 

0.13 



0.61 



+0.04 
-0.21 



t 



0.23 



+0.02 

0.05 



0.51«0.60) 44.2/53 0.60+°-°l 0.23+°-° 2 6 



0.31 



+0.03 



0.42 



+0.03 
0.18 



Notes: Spectral models: M=MEKAL thermal plasma (solar abundances) and PL=power-law continuum model. "Absorption column in 
units of 10 21 cm -2 . h Observed fluxes in the 0.3-10keV band, in units of 10 -13 erg s — 1 cm -2 . "Observed and unabsorbed luminosities in 
the 0.3-10 keV band, in units of 10 41 erg s 1 (assuming a distance of 56.7 Mpc). tSame hydrogen column as applied to the previous 
spectral component. Errors correspond to 90 per cent confidence limits for parameter of interest. The best-fit models are highlighted in bold. 
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Figure 2. Top left: XMM spectra and M+M+PL+GAUSS model for NGC 7771. PN data points and model fit are shown in black; those of the MOS data are 
shown in grey. Top right: The unfolded M+M+PL+GAUSS model and PN spectrum to illustrate the contributions from the model components. Bottom left: 
XMM spectra and M+M+PL model for NGC 7770 and bottom right: the unfolded PN spectrum. The model components are denoted by dashed (power-law), 
dotted (cool MEKAL) & dot-dashed (warm MEKAL) lines. 
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tectable in data of sufficient spatial and spectral quality for point 
source subtraction (e.g. NGC 3256. lLira et al JI2002T . we discount 
the M+M models as the hot plasma temperatures (kT ~ 19keV) 
are unrealistically high. Table|2|shows the results of the M+PL fits. 
We have included models where both absorbing components were 
were tied together (model 1) and where they were fitted indepen- 
dently (model 2), although there is little difference between them. 

However, although this simple model does adequately de- 
scribe the data, we have also tried fitting the spectra with a more 
complex three component model (M+M+PL), as it is unrealis- 
tic to assume that the thermal plasma will possess one distinct 
temperature. Indeed, these systems are expected to possess ther- 
mal pl asmas with a range o f temperatures spanning T ~ 10 J — 
10 s K jStrickland & StevenjEoOd) . This three-component model 
yields plasma temperatures of fcT=0.4/0.7 keV and a slightly softer 
power-law with T=1.7, with a 98.7 percent improvement over the 
M+PL model. We therefore adopt this as our best-fit model (model 
3), and Figure [2] shows the EPIC data plus residuals (top left) and 
the additive model components with the PN data (top right). This 
model implies that, while there is no additional absorption asso- 
ciated with the cool MEKAL component above the Galactic value 
(4.35x10 cm - 2 . lDickev & Lockmarll990ll . the power-law con- 
tinuum and warm thermal plasma component are much more heav- 
ily absorbed than in the simple M+PL models, although there are 
large uncertainties on these absorption parameters. Note that we 
were unable to improve the fit any further with the addition of more 
MEKAL components. 

Table|3|shows the relative contributions of the power-law and 
MEKAL components to the total X-ray emission in NGC 7771. The 
emission is dominated by the power-law component, indicating the 
presence of multiple XRBs or a low-luminosity AGN (LLAGN). 
The spectra for this source do show evidence of possible line emis- 
sion at 6-7 keV, which could be evidence for AGN activity (neutral 
Fe- K at 6.4 keV) or e mission from type Ib/IIa SNRs (Fe 6.7 keV, 
e.g. B ehar et alj200ll) in the central starburst region. To investigate 
this, we have attempted to fit a narrow gaussian line to the power- 
law continuum of the best-fit M+M+PL model. If all model param- 
eters were left free, we were unable to constrain the line energy. 
However, restricting the line energy to fit between 6 and 7 keV and 
freezing all other model components resulted in a line energy of 
6.24±g:fg keV with an equivalent width (EW) of 311 eV, with a 90 
per cent upper limit of 857 eV (plotted in Figure|2| top). If the line 
energy was frozen at 6.4 keV (i.e. at the expected position of a neu- 
tral Fe-K line from an AGN), the measured EW was 323 eV with 
an upper limit of 915 eV. The addition of this line only improves 
the fit statistics by ~ 1.7a; however its implications to the possi- 
bility of the presence of an LLAGN in NGC 7771 are discussed in 
section l5~2l 



4.1.3 NGC 7770 Spectral Analysis 

The X-ray spectra of the dwarf companion NGC 7770 were fit- 
ted with the same set of spectral models, and the results are also 
shown in Table [2] These data are well-fit with a M+PL model, al- 
though a statistically better fit is obtained by fitting the two absorb- 
ing columns separately (model 2). It is, however, unrealistic to have 
negligible absorption for the power-law component and substantial 
absorption for the thermal component as this model implies, since 
we assume that they originate from the same (star-forming) areas 
of the galaxy (i.e. the central and disc regions). Given the more 
limited statistics of the data, we therefore assume that both corn- 



Table 3. Fluxes, luminosities and percentage of total emission from the 
best-fit models of the central source in NGC 7771 and NGC 7770. 



Component 


Fx 


Per cent 


Lx 


Per cent 






of total 




of total 


NGC 7771 Central Source M+M+PL 


Cool MEKAL 


0.44 


16 


0.18 


8 


Warm MEKAL 


0.41 


15 


0.87 


39 


Power-law 


1.93 


69 


1.19 


53 


Total 


2.78 


100 


2.24 


100 


NGC 7770 M+M+PL 


Cool MEKAL 


0.06 


10 


0.10 


24 


Warm MEKAL 


0.22 


37 


0.14 


33 


Power-law 


0.32 


53 


0.18 


43 


Total 


0.60 


100 


0.42 


100 



Notes: Observed fluxes in the 0.3-10 keV band, in units of 10 

erg s — 1 cm -2 . Unabsorbed luminosities in the 0.3-10keV band, in 

units of 10 41 erg s — 1 (assuming a distance of 56.7 Mpc). 



ponents have the same absorbing column (model 1, fcT=0.5keV, 
T=2.4). 

However, as in the case of the central source in NGC 7771, 
the M+PL fit can be improved with the addition of a second MEKAL 
component. Again, given the more limited statistics for this galaxy, 
a model with separate absorbing columns for each component 
yields unrealistic values of Nh, and so we have applied one col- 
umn to all the components (model 3). Although partially uncon- 
strained, we are able to fit a low-temperature thermal component 
with kT=0. 15keV to the data, which improves the fit statistics at 
the 98.9 per cent (2.5<r) level. The data and model are plotted in 
Figure|2| and Table|3|shows the contributions of the power-law and 
MEKAL components to the total X-ray emission for the M+M+PL 
model. In contrast to the central source in NGC 7771, the thermal 
components contribute (in total) a larger fraction of the unabsorbed 
luminosity than the power-law component. This agrees with the 
soft nature of the galaxy in the XMM EPIC images (Figure^. Since 
this galaxy does have such a soft spectrum, we have attempted to fit 
it with thermal components only. A two-temperature plasma model 
(M+M) gives an unsatisfactory fit to the data (fcT=0.2/3.1 keV, 
xt ~ 1.4). If however we replace the soft power-law in model 3 
with a MEKAL component to fit the high energy end of the spectrum 
(i.e. M+M+M), this does yield a good fit (xt ~ 0.9) with plasma 
temperatures of kT — 0.16/0.56/4.00 keV. However this inter- 
pretation is problematic as the absence of a non-thermal component 
implies that there are no bright XRB systems present in NGC 7770. 

4.1.4 ULX Properties 

The two bright discrete sources located at either end of the galactic 
bar in NGC 7771 (X-l at a=23 h 51 m 22.3 s , S=+20° 06' 37" ; X-2 
at a=23 h 51 m 27.2 s , <5=+20° 06' 52" ) both have X-ray luminosi- 
ties in the ULX regime (> 10 39 erg s _1 ), and both have sufficient 
counts for simple spectral fitting (> 300 in the combined MOS & 
PN data). 

To begin with, we fitted the spectra of both sources with 
simple absorbed single-component spectral models: PL, MEKAL, 
DISKBB and BBODY. Of these, only the absorbed PL models gave 
statistically acceptable fits with \t < 2 for both sources (Ta- 
ble[4j. However, since these fits still gave a xt well above 1, and 
both showed evidence for a soft excess (Figure [3] left), we at- 
tempted to improve the fits with the addition of a soft component 
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Table 4. Spectral fitting results for the two ULXs in NGC 7771. 



N H kT/T in T x 2 /dof 
(keV) 


A x 2a 


l-P(F-test) 6 

(%) 




Obs 


d 

Unabs 


Flux Frac e 
Soft Hard 


NGC 7771 X-l 


Model 1: PL (WABS*PO) 

1.03±°;|^ 1 -72±o'.24 26.6/16 






5 11+ ' 72 
J11 -1.22 


| n 7 +0.28 
i,! "-0.47 


9 9S + 0.31 
z - z -0.65 






Model 2: PL+M (WABS*(PO+MEKAL)) 
2.36^;« 0.281°;^ L»±°;£ 14.2/14 


12.4 


98.8 


r c,+0.60 
-'■-' -2.51 


, 17 +0.23 
z - lz -0.97 


, qc+0.25 
z -"-2.23 


0.10 


0.90 


NGC 7771 X-2 


Model 1: PL (WABS*PO) 

<1.21 - l-89lo 25 21.7/13 






a n +0.66 
-'■ lz '-1.22 


1 20+ - 25 

lzu -0.47 








Model 2: PL+DISKBB (WABS*(PO+DISKBB)) 
3LM^ W6lg;g 15.9/11 


5.8 


82.0 


3.30t«;« 


1 27+ - 17 
1 - Z/ -0.99 




0.20 


0.80 


Model 3: PL+BBODY (WABS*(PO+BBODY)) 
3 ,3+4.04 1 2+0-04 1 7 3+o.3i , 60/11 

J - ZJ -1.61 u - lz -0.04 -0.54 ro.U/11 


5.7 


81.1 




, 9J-+0.18 
1,ZJ -1.04 


3.50^3 17 


0.18 


0.82 


Model 4: PL+M (WABS*(PO+MEKAL)) 
<1.98 <'.58to.3o L64 ia29 16/7/11 


5 


76.1 


3 3 4 +0.69 
■'•■"-1.30 


1 - zy -0.50 


, ,,+0.23 
1 -'-'-0.59 


0.08 


0.92 



Notes: Spectral models and parameters as in Tablel2lexcept: BBODY=blackbody model and DISKBB=multicolour disc black body 
model. "Improvement in the x 2 statistic over the single component model power-law (PL) fit, for two extra degrees of freedom. 
b F-test statistical probability of improvement of fit over single component model (PL). c Obser\>ed fluxes in the 0.3-10 keV band, in 
units of 10 — 14 erg s _1 cm -2 . d Observed and unabsorbed luminosities in the 0.3-l()keV band, in units of 10 40 erg s -1 (assuming a 
distance of 56.7 Mpc). Traction of total flux in the hard (PL) and soft model components over 0.3-10 keV band. The best-fit model 
for each source is highlighted in bold. 




Energy (keV) Energy (keV) 



Figure 3. XMM spectra of the two ULXs in NGC 7771. Top: PL (left) and PL+M (right) models for X-l. Bottom: PL (left) and PL+DISKBB (right) models 
for X-2. PN data points and model fit are shown in black; those of the MOS data are shown in grey. 
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NGC 7771 X-1 



NGC 7771 X-2 



10* 2xl0 4 

Time after start of observation (s) 



3xl0 4 



Figure 4. XMM-Newton light curves for the two ULXs in NGC 7771, ob- 
tained by summing the counts from the 3 EPIC cameras. The mean count 
rate is shown as a dashed line and error bars correspond to lc deviations 
assuming Gaussian statistics. 



(MEKAL/DIS KBB/B B ODY) , the results of which are also shown in 
Table|4] In the case of X-1, the best fit was achieved with a statisti- 
cal improvement at the 98.8 percent (2.5a) level with the addition 
of a soft MEKAL thermal component with kT ~ 0.3 keV and the 
resulting fit is shown in Figure[3](top right). Even though the addi- 
tion of either DISKBB & BBODY components also improved the fit, 
the resulting temperatures were very low (~ 0.07 keV), and both 
models had large absorption columns which resulted in unrealis- 
tically high unabsorbed luminosities (Lx > 10 42 ergs -1 ). For 
X-2, the PL fit was only improved with the addition of a soft com- 
ponent at the ~ 1.3a level (marginally best-fit with a PL+DISKBB 
model, Figure|5] bottom right), although this could simply be a con- 
sequence of the more limited photon statistics in this fainter source. 

Each source was tested for variability over the duration of 
the XMM-Newton observation using both x 2 an d Kolmogorov- 
Smirnov (K-S) tests. Firstly, the \ 2 test was use d to search for 
large amplitude variability in binned data against the hypothesis 
of a constant count rate, and for this purpose we derived short-term 
light curves for each source (Figure [4j- The data from the 3 EPIC 
cameras were co-added for improved signal-to-noise ratio, and tai- 
lored so that each bin had at least 20 counts after background sub- 
traction, resulting in temporal resolutions of 2000s for each source. 
With this test, X-2 showed some sign of short-term variability with 



Xv ~ 1-9, equivalent to a 97 percent (~ 2.2a) probability of vari- 
ability, while X-1 showed no evidence of variability (xt ~ 1.1). 
Secondly, the K-S test was used to search for smaller gradual small 
amplitude variations by comparing the observed background sub- 
tracted cumulative photon arrival distribution with the expected dis- 
tribution if the flux was constant. These tests were conducted using 
background subtracted PN light curves with a time-resolution of 
1 s. Here, it was X-1 that showed evidence of variability at the 97 
percent (~ 2.2a) level, while the data for X-2 showed none (< la). 



4.1.5 Comparison with previous observations 



using both 



NGC 111 HQ has previously been observed in X-rays using 
Einstein and ROSAT iFabbiano et alJll992t |Rej>haeli et alJ ll995t 
iDavies et ai]|l997l) . IBavies et alj Jl997l) estimated a total (0.16- 
3.5 keV) luminosity of Lx ~ 2.3 x 10 41 ergs" 1 for the com- 
bined emission of NGC 7771 and NGC 7770 from the Ein- 
stein data (corrected for the distance used in this paper), assuming a 
Raymond-Smith (RS) thermal plasma model with fcT=1.7keV. In 
the same study, ROSAT PSPC images showed strong X-ray emis- 



sion from the galaxy's centre, plus extended emission through- 
out the disc. Spectral analysis of the PSPC d ata was res tricted 
due to the low number of counts (~100), and lDavies et alj <1997l) 
were unable to distinguish between power-law and RS thermal 
plasma models. The power-law model gave an X-ray luminosity 
of 6.8 x 10 40 ergs -1 (0.1-2.4 keV) for the combined emission 
of NGC 7771 and NGC 7770, which is consistent with the total 
observed luminosity of the system in the same band in the XMM- 
Newton observation (Lx ~ 7.6 x 10 40 erg s _1 ). 



4.2 NGC 2342/2341 

4.2.1 X-ray morphology 

In Figure|5]we show soft (0.3-2 keV, top) and hard (2-10 keV, mid- 
dle) XMM-Newton EPIC images of the NGC 2342/1 pair of galax- 
ies. The images are the sum of all 3 detectors, and are convolved 
with a la (4 arcsecond) 2-D Gaussian mask and displayed with 
a logarithmic scale. The intensity minima/maxima correspond to 
surface brightnesses of 8.5 x 10~ 6 /5.4 x 10 -4 ct s -1 (soft) and 
8.5 x 10~ 6 /2.7 x 10~ 4 ct s -1 (hard). The contours correspond 
to 0.7,1.4 & 3.4 x 10~ 4 ct s" 1 (soft) and 0.3,0.7 & 1.3 x 10~ 4 
ct s _1 (hard). We do not resolve any detailed structure in either 
galaxy with XMM-Newton, but the companion galaxy NGC 2341 
does show an interesting elongated morphology in the hard X-ray 
image. Figure [5] (bottom) shows a DSS blue optical image of the 
system overlaid with broad-band (0.3-10 keV) X-ray contours cor- 
responding to surface brightnesses of 0.9,2.0 & 4.7x 10~ 4 ct s~\ 



4.2.2 NGC 2342 Spectral Analysis 

The integrated X-ray spectra for each of the galaxies were fitted 
using the same method described in section l4"Tl and the results are 
shown in Table [5] Although the spectra of NGC 2342 can be well- 
fit with a M+PL model, the addition of another MEKAL component 
improves the fit at the 98.3 percent (2.4a) level (model 3). We have 
had to use a single absorbing column for all components, as ap- 
plying separate columns to each component resulted in unrealistic 
values of N H (as with NGC 7770, see section l4"l~3l . This M+M+PL 
model (plotted in Figure© top) yields kT ~ 0.3/0.7 keV, T ~ 2.0, 
and an absorption colu mn of ap proximatel y twic e the Galactic 
value (8.93 x 10 20 cm~ 2 . |Pickev & Lockmarjjl990h resulting in an 
unabsorbed luminosity of Lx ~ 2.7 x 10 41 erg s _1 . Note that in 
this case, the two-temperature thermal plasma model (M+M) does 
not give a satisfactory fit to the data (xt ~ 1-5). The relative contri- 
butions of the spectral components for the best-fit model are shown 
in Table© the power-law component is very dominant in this case. 



4.2.3 NGC 2341 Spectral Analysis 

The spectra of this galaxy are best-fit with a simple M+PL model; 
in this case, the addition of another MEKAL component did not im- 
prove the fit above the 70 per cent level. As with NGC 7770, even 
though a superior fit is obtained with a M+PL with separate absorb- 
ing Nh columns for the two components (Table [5] model 2), the 
statistics for this source are too limited to separate the absorption 
columns in this way. Again, it is unrealistic to have a cool ther- 
mal plasma component which is more absorbed than the power- 
law component, and we therefore adopt the M+PL model with a 
single absorbing column as the best-fit model (model 1, plotted in 
Figure© bottom). The data are also well-fit (xt ~ 1.1) with the 
two-temperature thermal plasma (M+M) model with temperatures 
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Table 6. Fluxes, luminosities and percentage of total emission from the 
best-fit models of NGC 2342 and NGC 2341. 




Component 


Fx 


Per cent 


Lx 


Per cent 






of total 




of total 


NGC 2342 M+M+PL 


Cool MEKAL 


0.19 


6 


0.32 


12 


Warm MEKAL 


0.15 


5 


0.16 


6 


Power-law 


2.64 


89 


2.25 


82 


Total 


2.98 


100 


2.73 


100 


NGC 2341 M+PL 


MEKAL 


0.34 


23 


0.79 


44 


Power-law 


1.17 


77 


1.00 


56 


Total 


1.51 


100 


1.79 


100 



Figure 5. XMM-Newton EPIC soft (0.3-2 keV, top) & hard (2-10 keV, mid- 
dle) X-ray images of NGC 2342/1 with surface brightness contours over- 
laid. Bottom: DSS optical blue image of NGC 2342/1 with broad-band 
(0.3-10 keV) X-ray contours o verla id. See text for further details. 
© 2004 RAS, MNRAS 0Q0.HTH41 



Notes: Observed fluxes in the 0.3-10keV band, in units of 10 — 13 
erg s — 1 cm -2 . Unabsorbed luminosities in the 0.3-10keV band, in 
units of 10 41 erg s _1 (assuming a distance of 71 Mpc). 



of kT ~ 0.5/7.8 keV (not shown in Table|3), but we favour the 
power-law for the hard component, interpreted as emission from 
a population of XRBs. In this case, the contributions to the unab- 
sorbed luminosity from the two components in the best-fit model 
are comparable (Table |6j, with an enhanced percentage from the 
MEKAL component compared to NGC 2342. 



4.2.4 Comparison with previous observations 

The only other X-ray observati on of this system was performed 
with the ROSAT PSPC jHenriksen & CousineaiJll999h . However, 
the spectrum was extracted using a circular region of 3.4 arcminutes 
radius, which included the emission from both galaxies. These data 
were fit with RS thermal plasma models in the 0.25-2 keV range, 
giving a best-fit temperature of kT ~ 1.1 keV and a luminosity of 
Lx ~ 1.9 x 10 41 erg s _1 (corrected for the distance used in this 
paper). The total X-ray luminosity of both galaxies in the XMM data 
in the same energy range is Lx ~ 1-2 x 10 41 erg s _1 ; this ~ 35 
percent reduction is likely to be due to the large extraction region 
used in the analysis of the PSPC data, rather than a genuine change 
in the X-ray output of the galaxies themselves, although the latter 
possibility cannot be excluded. 



5 DISCUSSION 

5.1 The origin of the X-ray spectral components 

iPersic & RephaeliH2002l) have recently used the spectral properties 
of M82 and NGC 253 as observed with BeppoSAX together with 
stellar-population evolutionary models to quantify the X-ray spec- 
tral components of SBGs. Contributions to the X-ray flux in the 
0.3-10 keV band are expected to come from XRBs, SNRs, Comp- 
ton scattering of ambient FIR photons off supernova-accelerated 
relativistic electrons, diffuse thermal plasma and a compact nucleus 
in the form of a starburst or AGN. At low energies (<2 keV), the 
dominant component is expected to be diffuse thermal plasma re- 
sulting primarily from shock heating via the intera ction of the low- 
densi ty wind with the ambient high-density ISM I Strickla nd et alj 
2000). Above 2keV, they predict that the spectrum will be domi- 
nated by power-law emission from bright neutron star XRBs, plus 
a possible contribution from non-thermal Compton emission or an 
AGN if present JPersic & Rephaeiil20o3) . 
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Table 5. Spectral fitting results for NGC 2342/2341. 



PL 


Mi 


M 2 


X 2 /dof 


Fx" 


Lx 


c 


N H a r N H a 


kT(keV) 


N H a kT{keV) 






Obs 


Unabs 


NGC 2342 


Model 1: M+PL (WABS*(PO+MEKAL)) 

^t° f 6 2m+_° .\l t 


n ,o+0.06 
U - J -0.05 




144.7/158 


2 Q4+° 22 
-0.28 


, 77 +0.13 

L '"-o.vr 


9 go+0.13 
Z - 8z -0.52 


Model 2: M+PL (WABS*PO+WABS*MEKAL) 
1.70«-j 4 4 2.02^;^ <L46 0.62^™ 




144.1/157 


296 +o.io 

7 -0.33 


, 7 n+0.06 
1 - ' -0.20 


2 4 7 +0.03 
Z -^'-0.46 


Model 3: M+M+PL (WABS*(PO+MEKAL+MEKAL)) 

, or+0.41 1 qq+0.13 + n?q+ - 09 
1 - 83 -0.36 1 - yv -0.10 T "•""'-0.05 


t 74+ 19 

1 "-'^-0.13 


137.3/156 


2OS+0.23 


U0±°;« 


2 7^+0.17 



NGC 2341 



Model 1: M+PL (WABS*(PO+MEKAL)) 

2 83+ 1 66 182+ 014 t 30+ 09 - - 92 3/87 1 51+° 15 91+ 09 179+°' 
Model 2: M+PL (WABS*PO+WABS*MEKAL) 

, 07+1.71 ,70+0-14 S 62+ 2 - 01 029+ ' 05 - - 90 6/86 1 54+ - 11 n 93+ ' 06 201+°' 

l- s/ _0.84 1/u -0.18 J - DZ -1.60 uzy -0.07 yU.O/SO L -->^-0.35 uyJ -0.21 zul -l. 

Notes: Spectral models and parameters as in Tabled "Observed fluxes in the 0.3-10keV band, in units of 10~ 13 erg s 1 cm . 
h Observed and unabsorbed luminosities in the 0.3-10 keV band, in units of 10 41 erg s — 1 (assuming a distance of 71 Mpc). tSame 
hydrogen column as applied to the previous spectral component. Errors correspond to 90 per cent confidence limits for parameter of 
interest. The best-fit models are highlighted in bold. 




Figure 6. Top left: XMM spectra and M+M+PL model for NGC 2342. PN data points and model fit are shown in black; those of the MOS data are shown in 
grey. Top right: The unfolded model and PN spectrum to illustrate the contributions from the model components. Bottom left: XMM spectra and M+PL model 
for NGC 2341, and bottom right: the unfolded PN spectrum. The model components are denoted by dashed (power-law), dotted (cool MEKAL) & dot-dashed 
(warm MEKAL) lines. 
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These predictions have in part been confirmed by XMM- 
Newton and Ch andra obs ervations of l o cal S BGs, e.g. M82 
jMatsumoto et alJ l200lb IStevens et alJ l2003h. NGC 253 
JPietsch et alJ 1200 it IStrickland et al J 1200(1 IStrickland et alj 
20021) NGC 3256 & NGC 3310 <Lira et al J 120021: 1 Jenkins etafl 
2004b) and NGC 4038/4039 (the Antennae: e.g. iFabbiano et alJ 
2001). Their broad band spectra can indeed be fit with single- 
or multi-temperature thermal plasmas with temperatures ranging 
between ~0.2-0.9keV, plus a hard power-law tail at energies 
>2keV. The interpretation of the origin of the power-law com- 
ponent must however be revised; the high spatial resolution of 
Chandra has shown that the hard emission can be attributed 
to a few dominant luminous compact sources which are most 
likely XRBs with black-hole rather than neutron star primaries 
(based on their high luminosities). Spectacul ar examples of th e 
detection of such sources include the Antennae iZezas et al. 2002), 
the Car twheel jGaoetal J 120031: IWolter & Trinchieril 120041) and 
Arp299 fcezas et all2003l) . 

Our spectral results for NGC 7771/0 and NGC 2342/1 are con- 
sistent with the other observations of nearby SBGs listed above. For 
both NGC 7771 and NGC 7770, we are able fit the thermal spec- 
tra with a two-temperature MEKAL models with fcT=0.4/0.7 keV 
and A;T=0.2/0.5 keV. In the case of NGC 7771, the trend of increas- 
ing absorption with increasing plasma temperature is similar to that 
found in other SBGs, and is evidence that the hotter component is 
likely to originate from supernova-heated gas in the most obscured 
central disc and starburst regions. A clear illustration of such a sce- 
nario is found in NGC 3256, where the hotter of the two plasmas 
can be spatially associated with the central ~ 2 kpc absorbed re- 
gion of the galaxy, with the cooler component dominating in the 
outer regio ns most lik ely in the form of an out-flowing galactic 
wind iLira et al.fl2 002). The complexity of the thermal models we 
are able to fit is highly dependent on the quality of the data, and we 
must bear in mind that these models are an over-simplification of 
the gas profiles in these galaxies i.e. it is physically unrealistic for 
there to be such distinct thermal plasma temperatures. This com- 
parison of spectral fits of differing quality also highlights another 
point; with good enough statistics we can separate out the absorb- 
ing columns associated with each thermal component, whereas this 
type of modelling can lead to unrealistic values of Nh when fit to 
spectra with insufficient counts (i.e. NGC 7770/2342/2341). 

In addition to the soft thermal components, we can use the 
information gleaned from Chandra observations of other SBGs to 
interpret the origin of the hard (> 2keV) integrated emission we 
detect from NGC 7770, NGC 2342 and NGC 2341, which are 
spatially unresolved by XMM-Newton. In these galaxies, there is 
no evidence for hidden luminous AGN activity, which can be de- 
tected at X-ray energies by the presence of an absorbed power-law 
(Nh > 10 22 cm -2 ) and/or a strong neutral iron line at ~6.4keV 
from the nucleus. In these cases, therefore, the power-law emission 
is almost certainly coming from populations of XRBs. However, 
there is some evidence for AGN activity in NGC 7771 (see sec- 
tionl5~2l. 

5.1.1 The ULXs in NGC 7771 

The two bright extra-nuclear X-ray sources in NGC 7771 are 
co-spatial with the ends of the galactic bar (see section 14.1.11 . 
This is consistent with models that show that gas can accumulate 
at the tips of b ars due to the corotation o f the bar structure 
with the disc lEnglmaier & Gerhard Il997l) . thereby creating 
star-forming regions where XRB sources can form. A few 



other examples of X-ray sources at bar-ends have also been 
observed, for example in ROSAT obse rvations of t h e barred s pirals 
NGC 1672 iBrandt. Halpern. & Iwasaw al Il996t Ide Narav et"ail 
l200d) and ROSAT/Chandra observations of NGC 4 303 
iTschoke. Hensler. & Junkesi2000l : ljimenez-Bail6n et all2003l) . 

The X-ray spectral fits to the ULXs are in the form of a power- 
law continuum plus a soft excess. X-l is well fit with a PL+M 
model (r ~ 1.6, kT ~ 0.3 keV, xl ~ 1-0), while the power- 
law fit to X-2 is marginally improved (although only at the 82 per 
cent level) with a DISKBB component (T ~ 1.7, T in ~ 0.2 keV, 
xl ~ 1-4), although normal BBODY and MEKAL components also 
give improved fits over a power-law only model. Th ese types of 
spectral fits are commonly seen in bright ULX s (e.g. Mill er et alJ 
l2003t I Jenkins et aljExM^t iRoberts et alJl2004l) . The physical na- 
ture of ULXs has been the topic of much debate in recent years, 
but XMM-Newton and Chandra observations are greatly improving 
our understanding of this phenomenon. The observed link between 
ULXs and episodes of high star formation, coupled with their spec- 
tral and timing properties suggest that the majority are likely to be 
high-mass XRBs going through a period of thermal-timescale mass 
transfer l Kins 2004). However, alternate explanations invoke, for 
examp le, the presence of i ntermediat e-mass black h oles (IMBHs, 
e.g. IColbert & Mushotzkvlll999t iMiller et alJl2003l) . Indeed, the 
possible presence of a such a soft accretion disc component in X-2 
is consistent with the IMBH scenario, based on the argument that 
black hole mass scales inversely with the accretion disc temperature 
(Ti n oc Af" 1//4 , e.g. Makishima et al. 2000). However, the data are 
photon-limited leading to ambiguous results (i.e. that the soft emis- 
sion can also be well-modelled with normal blackbody and thermal 
plasma models), and the DISKBB component is not is not signifi- 
cant enough for us to claim this as an IMBH candidate at present. 
The soft component in the spectrum of X-l, again at low spectral 
resolution, can be modelled by a thermal plasma which is likely 
to originate in the star-forming region where the source resides in 
a similar manner to the thermal plasma component in the central 
source in NGC 7771. The short-term variability detected in both 
sources (section l4.1.4> strongly suggests that their X-ray emission 
is dominated by single accreting XRBs, although high spatial reso- 
lution imaging with Chandra would be required to confirm this. 

5.2 Does NGC 7771 host a low-luminosity AGN? 

NGC 7771 possesses a large-scale bar and circumnuclear ring, 
both of which are likely to be a consequence of disc instabili- 
ties b rought on by its interaction with NGC 7770 ISmifh et alJ 
1999). The presence of a bar plays an important role in a 
galaxy's evolution; bars can extract angular momentum from the 
ISM in the disc through gravitation torques, causing gas to be 
driven into the inner few hundred pc of the galaxy, thereby 
fuelling a burst of star formation and possibly an AGN (e.g. 
iBarnes & Hernquistlll99rl) . Observational evidence clearly shows 
that nuclear starburst activity is enhanced in galaxies with bars, 
often in the form of circumnuclear rings resulting from gas ac- 
cumulating in the vicinity of in ner Lind blad reso nances (ILRs; 
iKnapen. Perez-Ramirez. & Lainei2002Hjogeel2004l and references 
therein). However, there is conflicting evidence with regards to 
the incidence of AGN in barred galaxies; while some studies 
find a higher fraction of bars in S eyfert galaxies compared with 
'normal' (non-active) galaxies (e.g.lKnapen. Shlosman. & PeletieJ 
2000; Laurikainen, Salo, & Buta 2004), other results show that 
bars have n egligible effe ct on the presence of an AGN (e.g. 
iMulchaev & Reganll997tlHo. Filippenko. & Sargentlll997iJ) . This 
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may be explained by the fact that models of gas flow in barred 
galaxies show that gas driven toward the nucleus can be interrupted 
by a nuclear ring, which may prevent the fueling of an AGN, al- 
though sufficient stellar material may be present in the immediate 
vicinity of galactic nuclei to fuel LLAGN iHo et alll997l and ref- 
erences therein). However, since the duty cycles of AGN are ex- 
pected to be a factor 10-100 times shorter than the lifetime of bars 
(~ 1 Gyr), we do not expect all barred galaxies to show signs of 
AGN activity ljpgeel2004l) . 

The lack of a correlation between bars and AGN may 
in part be explained by obsc uration, i.e. 'hidden' AGN. 
iMaiolino. Risaliti. & Salvatil 1 19991) reported that highly obscured 
AGN (e.g. Compton-thick Seyfert 2; N H > 10 24 cm~ 2 ) are 
found preferentially within strongly barred galaxies, whereas less 
obscured AGN are located in hosts with weak bars or no bar 
present. This correlation indicates that gas and dust driven to- 
wards the nucleus creates an obscuring screen which prevents 
us from viewing the nucleus directly. In addition, starbursts are 
naturally expected to accompany the formation of an AGN due 
to the large amounts of gas in-flowing into the central regions, 
and it is feasible that the starburst itself (if extended) can ob- 
scure the AGN, resul ting in optical spectra that appear to us as 
HII -type objects iFabian e t al,||l99 8f). Clear exam ples of this phe- 
nomena are NGC 494 5 iGuainazzi etail EoOQ^ and NGC 3690 
iDella Ceca et al. 2002), both of which show evidence of AGN ac- 
tivity at X-ray wavelengths, although their optical and IR spectra 
are dominated by the starburst component and are classified as HII. 

IMaiolino et alJl20ol report that many hidden AGN are found 
in galaxies with 10 11 < Lir/Lq < 10 12 , the range in which 
NGC 7771 lies. Even though no spectral signatures of AGN activ- 
ity have been found at other wavelengths, a comp act nuclear s ource 
is detected at radio and IR wavelengths ISmith et alJll999^) . The 
presence in the new X-ray data of both a hard absorbed power- 
law com ponent with a slope of F ~ 1.7 typically seen in Seyfert 
galaxies (Nandra & Pounds 1 19941) . plus a possible iron line at ~ 
6keV in the central source in NGC 7771 suggests that an ob- 
scured AGN may be present. We can test this hypothesis in sev- 
eral different ways. The absorption we measure in the X-ray data of 
Nh ~ 10 22 cm -2 is equival ent to Ay ~ 5.5 (assuming a sta ndard 
Galactic value of E B -v /N H JBohlin. Savage. & Drakell978l) . The 
optical colour excess measured by lVeilleuxetalll995l) implies an 
ex tinction consistent with this of Ay ~ 6.5. Using the relation 
of lElvis. Soltan. & Keell 1 19841) between the hard (2-10 keV) X- 
ray luminosity and broad Ha emission-line luminosity in LLAGN 
(Lx I 'L(Ha) ~ 40), we can determine whether we would expect 
to detect a broad Ha line from the nucleus if one were present, or 
if it would be entirely obscured. The 2-10 keV flux of NGC 7771 
of ~ 1.8 x 10~ 13 erg cm -2 s _1 predicts a broad Ha component 
with F[Ha) ~ 4.5 x 10" 15 erg cm~ 2 s _1 , which is a factor ~ 
6 less than the flux of the narr ow line detected in the most ob- 
scured region of the nucleus bv lDavies et all ll997l) . With 6 mag- 
nitudes of extinction, this flux would be reduced by a factor ~ 
250, so a broad Ha component would be beyond detection lim- 
its. There is also a UV source in the nucleus detected by IUE with 
a flux of ~ 2.5 x 10~ 12 erg cm -2 s _1 (see figure 8. lDavies et all 
1 19971) . However, this emission is likely to be coming from the star- 
burst ring (~ 6 ar csec diameter), where the obscuration is patchy 
ISmithet all 19991) . since the ratio of IR to UV flux on the spectral 
energy distribution (SED) of a factor of ~ 5 implies a maximum 
Ay of ~ 1 magnitude in the UV if the intrinsic SED were that of 
an AGN. This is inconsistent with the extinction measured in the 
optical in the central ~ 2 arcseconds of the galaxy. 



The final tentative piece of evidence for a LLAGN is the low- 
significance ( 1 J a) detection of a narrow Fe-Ka line at ~6 keV. We 
can predict the EW of such a lin e from an active nucleus using the 
empirical relation of IPage et alj 120041) . This describes the X-ray 
'Baldwin effect', in which the EW of the line decreases with in- 
creasing luminosity (EWoc ^(^lOfceV))' P 0SSlrJ ly resulting from a 
decrease in the covering factor of the a molecular torus around the 
AGN. The 2-10 keV luminosity of L x ~ 7 x 10 40 erg s _1 cor- 
responds to an Fe-Ka EW of ~ 240 eV, which is consistent with 
the measured (although partially unconstrained) EW of the line in 
the spectrum of the central source in NGC 7771. Another impor- 
tant point is that while the 90 per cent errors on the line energy 
(6.24^q' 2 q keV) are consistent with 6.4 keV emission, they do rule 
out the possibility that it originates from ~ 6.7 keV helium-like 
iron from the starburst component. Future X-ray observations with 
high signal-to-noise data at ~ 6 keV will be required to confirm the 
presence of this line, which could firmly establish the existence of 
a LLAGN in NGC 7771. 



5.3 The consequences of interaction 

We know that interactions between galaxies increase their lev- 
els of star formation as seen in the optical/IR (see section 0. 
To investigate how the X-ray properties of our interacting sys- 
tems relate to those of other galaxies, we have compared them 
with three different samples; isolated galaxies, interacting galax- 
ies and merger systems. We note that these are in no sense com- 
plete samples, but are simply representative samples for com- 
parison purposes. The X-ray data for the isolated galaxy sam- 
ple (M101, NGC 278, M74, NGC 1291, NGC 2681, NGC 3184, 
NGC 4314, M94, M83 & IC 5332) and the interacting sys- 
tem M51A/B ar e taken from t he Ch andra survey of nearby spi- 
ral galaxies of iKilgard et al] 120041) . The remainder are taken 
from XMM-Newton, Chandra or ROSAT IASCA observations in 
the literature. In addition to the galaxies st udied in this paper, 
the interacting systems are: The Cartwheel I Wolter & Trinchieri 
|2004 [); M82/M81/NGC 3077 IStevens et all l2003t ISwartz et al' 
2003[ lOtt. Martin & Waited 120031) . NGC 4485/90 iRoberts et al' 
2002), while the sample of merging galaxies comprises NGC 3690 
(Arp 299; IZezas. Georeant opoulos. & Ward Il998l). NGC 3256 . 
NGC 331 & NGC 4438/39 (tak en from Ijenkins et all l2004bl) . 
Aid 220 iMcDowell et all l2003l) . NGC 3921 & NGC 7252 
iNolan et al]l2004l) and NGC 4676 lRead E0031. For consistency 
between the datasets, these X-ray luminosities have all been con- 
verted to the 0.3-10 keV band. 

The FIR luminosities for each galaxy are calculated using 
FIR fl uxes d e rived f rom the 60 and 100 /xm flux densiti es of 
ISanders et ail 120031) . ISurace. Sanders. & Mazzarellal 120041) and 
the IRAS point s ource catalogue v2.0 (via the NED archive), us- 
ing the relation of Hel ou. Soifer. & Rowan-Robinsonl 1 1985): 

FIR = 1.26 x 10 -11 (2.58S6CV + S100J er S cm ~ 2 s_1 

The optical B-band luminosities are calculated using the ex- 
pression of iTullvl 19881) relating the blue apparent magnitude Bt 
and the distance D (in Mpc): 

logL B (L @ ) = 12.192 - 0AB T + 2logD 

The blue apparent magnitudes are mainly taken from the RC3 
catalogue Ide Vaucouleurs 1991), with the value for the Cartwheel 
taken from the ESO/Uppsala survey of the ESO(B) Atlas iLaubertsl 
1982). The radio luminosities are derived from the 1.4 GHz flux 
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Figure 7. Plots of the 0.3-10 keV X-ray luminosity against FIR (top), blue 
(middle) and radio (bottom) luminosities for the three samples of galax- 
ies. Isolated galaxies are marked with open squares, merging systems are 
marked with open circles, interacting galaxies are marked with filled trian- 
gles (with NGC 7771/0 and NGC 2342/1 marked with black triangles). See 
text for details of data sources. Note that NGC 7770 is not shown in the 
Lx/Lpm plot as there are no IR fluxes catalogued for this galaxy. 



densities of ICondon et alj ll998l) where available, whil e the re- 
maind er are derived from the 1.4 9 GHz flux densities from lCondonl 
Jl987l) and lCondon et alj <199fJ) . 

The 0.3-10 keV X-ray luminosities are plotted against these 
quantities in Figure Q for the three galaxy samples. The major- 
ity of the sample show a strong correlation between the X-ray lu- 
minosity and their multiwavelength properties. This is to be ex- 
pected, since the optical blue flux is consistent with the pres- 
ence of young stellar populations, and both FIR and radio lumi- 
nosities are strong tracers of the star formation activity associ- 
ated with those populations. The FIR emission from SBGs comes 
primarily from dust heating by UV flux from massive OB stars 
in the starburst regions, which is re-radiated in the IR, while 
the nonthermal 1.4 GHz radio flux is believed to arise from syn- 
chrotron radiation created by the acceleration of cosmic-ray elec- 
trons by supernova explosion s (plus a small fraction from SNRs 
thems elves) (e.g. iHelou et"allll985| : ICondon. Anderson. & Helot] 
1 1 99 if) . X-ray emission is also an excellent tracer of star forma- 
tion, since it arises from phenomena related to the end points of 
stellar evolution, mainly in the form of XRBs, SNRs and diffuse 
thermal plasma. A correlation between the FIR and ~0.2-4.5 keV 
X-ray emissi on in spiral an d bright I R galaxies was first seen us- 
ing E instein fDavid. Jo nes, & Forman||l992|: |Fabbiano & Shaplev 
20021). More recently. iRanalli. Comastri. & Settil J2003I) and 
Gilfanov, Grimm. & Sunvaev (2004) have demonstrated that there 
is a tight linear correlation in SBGs between their FIR, radio and 
soft (0.5-2 keV) and hard (2-10 keV) X-ray luminosities, which we 
extend here to include isolated (non-starburst) galaxies. When the 
co rrelation between th e FIR and 1 .4 GHz luminosities was reported 
bv lHelou et all £985), it was also noted that it held over a wide 
range of objects, from quiescent spirals to galaxies dominated by 
nuclear starbursts. This agrees with extension of the isolated spiral 
galaxies in the X-ray correlations shown here. 

These relations are continuous between the isolated, interact- 
ing and merging systems. However, while there is a clear segre- 
gation between the isolated and merging galaxies, the interacting 
sample span almost the entire range in luminosities at all wave- 
lengths, even extending up to the level of some of the merging 
systems. The three interacting galaxies that appear to deviate from 
these correlations with X-ray luminosities roughly an order of mag- 
nitude higher than the rest of the sample in the Lx-Lfir plot are 
M81, whose X-ray luminosity is dominated by its bright nucleus 
( Swart zet alJl2003T) . the Cartwheel, which has numerous bright 
ULXs in a star burst ring away f rom the obscured centre of the 
galaxy iWolter & TrinchierilEool) . and NGC 4485, whose X-ray 
output is dominated by one bright ULX iRoberts et all2002h . 

Ranalli et al. 1 2003) find that while some galaxies with known 
LLAGN have Lx-Lfir and Lx-L\aghz ratios exceeding those 
of SBGs, some do follow the starburst relation, demonstrating that 
their outputs at these wavelengths are dominated by the star forma- 
tion activity rather than the AGN component. An example of this is 
seen in NGC 3690 (Arp 299), where the emission from the buried 
AGN only contribu tes a very small fraction of the bolometric emis- 
sion of the galaxy jDella Ceca et all2002tlZezas et alJl2003h . It is 
clear from these plots that the galaxies studied in this paper (marked 
with black triangles) are among the brightest interacting galaxies in 
X-rays, and are at the high-luminosity end in all other properties, 
overlapping with some merging systems. In addition, these galax- 
ies closely follow the SBG relation, demonstrating that their X-ray 
and FIR luminosities are dominated by the star formation activity, 
including NGC 7771 where there is a question as to the presence of 
a LLAGN. 
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6 SUMMARY & CONCLUSIONS 

In this paper, we have studied the X-ray properties of two pairs 
of interacting galaxies, NGC 7771/0 and NGC 2342/1, in the 0.3- 
lOkeV energy band. For the first time, with the high quality spectra 
obtained with XMM-Newton we have detected the classic X-ray sig- 
natures of starburst activity in all four galaxies, in the form of soft 
multi-temperature thermal emission dominating below 2keV and 
hard power-law components representing the emission at higher en- 
ergies (2-10 keV). 

In NGC 7770, NGC 2342 and NGC 2341, we do not re- 
solve any X-ray structure, but their integrated emission can be 
well-modelled with multi-temperature thermal plasmas with tem- 
peratures range between kT ~ 0.2-0.7 keV, with the warm com- 
ponent possibly representing supemova-heated gas within the ab- 
sorbed discs and central regions and the cooler components possi- 
bly representing out-flowing material from the starburst interacting 
with the ISM. The power-law components have slopes of V ^4.8- 
2.3, and are likely to represent the integrated emission from pop- 
ulations of XRBs associated with the starbursts, similar to those 
seen in other nearby merger systems. In NGC 7771, we resolve 
the X-ray emission into a bright central source plus two ULXs lo- 
cated at the either end of its bar at the corotation radius. The soft 
X-ray spectrum of the central source is well-modelled by a two- 
temperature thermal plasma with kT ~0.4/0.7 keV, showing in- 
creasing absorption with increasing plasma temperature. However, 
the power-law component is relatively flat and absorbed (F ~ 1.7, 
Nh ~ 10 22 cm -2 ). This, coupled with a low-significance detec- 
tion of an emission line with an EW ~300 eV at ~6keV, is the first 
indication that NGC 7771 may host a LLAGN, perhaps fuelled by 
the inflow of material along the galactic bar. The two bright bar 
sources are well-fit with power-law plus soft excess models com- 
monly seen in ULXs, and both are variable above the 2a level, im- 
plying that their emission is dominated by single accreting XRBs. 

A comparison with other isolated, interacting and merging 
systems from the literature shows that the galaxies studied in this 
paper are among the brightest interacting systems at X-ray, FIR, 
B-band and radio luminosities, and are even comparable to some 
late-stage merging systems. They all follow the established lin- 
ear correlations for starburst systems between Lx-Lfir, Lx- 
Lb and Lx-Liaghz, demonstrating that their X-ray output is 
dominated by emission from the starburst components, including 
NGC 7771 which may host a LLAGN. These types of in-depth 
studies of interacting and merging systems at X-ray energies with 
XMM-Newton and Chandra can assist in our understanding of the 
evolution of galaxies as a whole, revealing whether such interac- 
tions increase the activity levels in galaxies in the form of pure 
starbursts, and whether interactions play a major role in supplying 
sufficient fuel to trigger AGN activity that may be obscured at other 
wavelengths. 
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